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Abstract 


The  equivalent  beam  model  is  widely  used  for  predicting  strength  and  vibration 
of  a  ship  hull  in  a  preliminary  analysis.  It  can  also  be  used  for  checking  results  in  a  large 
finite  element  model  and  for  parametric  studies  of  ship  behaviour.  This  method  treats  a 
ship  hull  as  a  series  of  prismatic  segments  connected  together.  Each  segment  has  its 
sectional  properties,  real  and  virtual  masses.  The  program  SCRAP  was  developed  at 
Defence  Research  Establishment  Atlantic  for  calculation  of  cross  sectional  constants  and 
estimation  of  mass  properties.  It  prepares  input  data  files  for  the  finite  element  programs 
VAST  and  TORSON  and  interprets  the  analytical  results.  SCRAP  can  be  only  used  for 
some  specific  sections  at  the  present  time.  For  arbitrarily  oriented  sections  it  may  give  an 
incorrect  shear  centre  and  warping  constant  and  thus  the  wrong  stress  distributions. 

This  report  presents  the  mathematical  derivations  of  the  equations  used  for  the 
calculation  of  cross  sectional  constants  and  stress  distributions  of  thin-walled  sections.  As 
an  improvement  over  the  current  SCRAP  program,  these  equations  are  applicable  to  any 
shape  of  cross  section,  both  open  and  closed,  and  are  independent  of  the  orientation  of 
the  cross  section.  A  computer-oriented  step-by-step  procedure  based  on  these  equations 
is  outlined.  Several  examples  are  also  presented  to  verify  the  procedure. 

Resume 

Dans  une  analyse  preliminaire,  le  modele  de  poutre  equivalente  est  largement 
utilise  pour  la  prevision  de  la  resistance  et  des  vibrations  d’une  coque  de  navire.  II  peut 
aussi  etre  utilise  pour  verifier  les  resultats  dans  un  modele  a  elements  finis  de  grande 
dimension  et  dans  des  etudes  parametriques  sur  la  tenue  des  navires.  Dans  la  presente 
methode,  on  considere  qu’une  coque  de  navire  est  constitute  d’une  serie  de  segments 
prismatiques  relies  ensemble.  Chaque  segment  possede  ses  proprietes  de  section,  sa  masse 
reelle  et  sa  masse  virtuelle.  Le  programme  SCRAP  a  ete  ecrit  au  Centre  de  recherches 
pour  la  defense  (Atlantique)  en  vue  du  calcul  des  constantes  de  section  transversale  et  de 
l’estimation  des  proprietes  de  masse.  II  prepare  les  fichiers  d’entree  pour  les  programmes 
a  elements  finis  VAST  et  TORSON  et  il  interprete  les  resultats  analytiques.  SCRAP  ne 
peut  pour  le  moment  etre  utilise  que  pour  certaines  sections  particulieres.  Dans  le  cas  de 
sections  d’orientation  arbitraire,  il  peut  donner  une  constante  de  centre  de  cisaillement  et 
de  gauchissement  erronee  et,  par  consequent,  des  distributions  de  contraintes  erronees. 

Ce  rapport  montre  comment  on  a  deduit  les  equations  mathematiques  utilisees  pour 
calculer  les  constantes  de  section  transversale  et  £tablir  les  distributions  de  contraintes  des 
sections  ^  pario  mince.  Ces  Equations  constituent  un  perfectionnement  par  rappon  an 
programme  SCRAP  actuel  :  elles  peuvent  etre  appliquees  a  toutes  les  formes  de  section 
transversale,  ouvertes  et  fermees,  et  elles  sont  inddpendantes  de  rorientation  de  la  section 
transeversale.  On  decrit  une  methode  pas  a  pas  m6canisee  basee  sur  ces  equations. 
Plusieure  exemples  permettant  de  verifier  la  methode  sont  aussi  presentes. 
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area  of  closed  section 

cross  section  area 
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1.  Introduction 


The  equivalent  beam  model  is  a  simplified  method  widely  us^i  for  predicting 
strength  and  vibration  of  a  ship  hull  in  preliminary  design.  It  can  Uso  be  used  for 
parametric  studies  and  for  checking  results  in  a  large  finite  element  model.  The 
equivalent  beam  method  treats  a  ship  hull  as  a  series  of  prismatic  segments  connected 
together.  Each  segment  has  its  own  sectional  properties,  real  and  virtual  masses.  By 
utilizing  different  types  of  beam  elements  and  numerical  methods,  stress  distributions  due 
to  various  loadings  and  natural  frequencies  of  the  structures  can  be  obtained. 

The  various  beam  elements  in  the  finite  element  method  are  derived  from  different 
assumptions,  with  or  without  shear  deformation  and/or  longitudinal  warping.  For  a  ship 
with  large  openings,  such  as  a  container  ship,  the  lowest  natural  frequency  of  coupled 
horizontal-torsional  modes  may  be  close  to  the  lowest  natural  frequency  of  flexural 
vibration  modes,  and  the  horizontal  propeller  forces  may  generate  large  torsional  moment. 
The  shear  deformation  of  the  cross  section  of  a  ship  hull  also  can  be  significant 
depending  on  the  depth  to  length  ratio  of  the  ship.  Thus,  the  required  sectional  properties 
for  the  subsequent  numerical  analysis  are  not  only  the  moments  of  inertia  and  centroid, 
but  also  the  shear  centre  and  torsion  and  warping  constants. 

The  program  SCRAP  was  developed  by  the  Structural  Mechanics  Group  of 
Defence  Research  Establishment  Atlantic  (DREA)  for  calculating  different  sectional 
constants,  mass  properties,  preparing  input  data  files  for  the  finite  element  programs 
VAST  fl]  and  TORSON  [2]  and  interpreting  the  analytical  results.  The  required  sectional 
properties  for  the  general  beam  element  in  the  finite  element  program  VAST  are  the 
moments  of  inertia,  the  torsion  constant  and  the  shear  centre.  An  additional  sectional 
constant  needed  for  the  program  TORSON  is  the  warping  constant.  SCRAP  is 
documented  in  References  3  to  6. 

The  program  SCRAP  is  efficient  and  user  friendly  but  can  only  be  used  for 
symmetric  sections  at  the  present  rime.  It  may  give  the  incorrect  shear  centre  and  warping 
constant  for  an  arbitrarily  oriented  section  and,  consequently,  invalid  results  for  the 
equivalent  beam  analyses.  Because  this  program  is  unable  to  provide  the  correct  shear 
centre  for  arbitrarily  oriented  sections,  the  distribution  of  transverse  shear  flows  is  also 
incorrect. 

The  available  literature  on  the  subject  of  calculating  sectional  constants  and  stress 
distributions  of  arbitrarily  oriented  thin-walled  sections  tends  to  be  incomplete  and 
problem  specific  [7,8],  This  report  presents  complete  mathematical  derivations  of  the 
equations  used  for  the  calculation  of  the  cross  sectional  constants  and  stress  distributions 
of  thin-walled  sections.  These  equations  can  be  used  for  any  shape  of  cross  section,  both 
open  and  closed,  and  are  independent  of  the  orientation  of  the  cross  section.  As  well,  a 
computer  oriented  step-by-step  procedure  based  on  these  equations  is  outlined.  Several 
examples  are  also  given  in  the  final  section  of  this  report  to  verify  this  procedure. 
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2.  Derivations 


Geometric  discontinuities  exist  at  large  hatch  openings  of  ship  hulls  and  on  the 
locations  where  open  and  closed  sections  mix.  They  create  relative  restraints  and  cause 
displacement  incompatibilities  between  segments  under  applied  torsion  loading,  thus 
generating  secondary  stresses  in  addition  to  the  primary  shear  stresses  and  altering  the 
overall  stress  distribution.  This  section  presents  the  mathematical  derivations  of  torsion 
induced  displacements  and  stresses  of  thin-walled  sections. 

Several  essential  assumptions  made  herein  for  all  of  the  derivations  follow: 

1.  The  beam  segment  is  a  prismatic  thin- walled  section.  The  term  "thin-walled  section" 
indicates  that  the  thickness  of  the  wall  Y  is  small  in  comparison  with  the  total  width  of 
the  cross  section  but  it  has  sufficient  thickness  so  that  local  buckling  is  not  a  problem. 
Thus,  the  shear  stresses  can  be  assumed  constant  through  the  thickness  of  the  wall  but 
may  vary  along  the  cross  section. 

2.  There  is  no  transverse  deformation  occurring  under  applied  loads.  The  shape  of  the 
cross  section  always  remains  unchanged. 

3.  The  shear  deformation  of  the  cross  section  is  caused  by  the  primary  (St-Venant)  shear 
stresses  only.  The  additional  shear  deformation  caused  by  the  secondary  (warping)  shear 
stresses  can  be  neglected.  For  an  open  thin  walled  section,  the  net  St-Venant  shear 
stresses  are  equal  to  zero  through  the  wall  thickness;  consequently,  there  is  no  shear 
deformation. 

The  sign  conventions  and  coordinate  systems  of  a  thin-walled  cross  section 
subjected  to  bending  and  torsion  for  the  following  theoretical  development  are  illustrated 
in  Fig.  1.  Two  coordinate  systems  are  defined;  a  Cartesian  coordinate  system  which  has 
x,  y  and  z  axes  through  the  centre  of  gravity  of  the  cross  section  and  a  curvilinear 
coordinate  system,  s,  which  coincides  with  centreline  of  the  wall  of  the  cross  section  and 
is  bounded  with  boundaries  s=0  and  s=b.  The  Cartesian  coordinate  system  has 
displacements  t)  and  w  along  the  x,  y  and  z  axes  and  a  rotation  (p  around  the  z  axis. 
The  curvilinear  coordinate  system,  on  the  other  hand,  has  a  displacement  \)  along  the  s 
axis.  The  positive  directions  of  forces  and  displacements  are  indicated  by  arrow  heads. 
One  should  keep  in  mind  that  the  sign  conventions  used  here  are  different  from 
Timoshenko’s  definitions  [9]. 

2.1  General  Considerations 

A  small  wall  element  cut  from  the  thin  wall  beam  as  shown  in  Fig.  1  has  an  angle 
a  between  the  s  and  x  axes.  The  displacements  between  two  coordinate  systems  can  be 
related  to  each  other  through  geometric  orientation  as: 
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Compatibility: 

3w  3o 

ir*ir ' r 


„  -N 


di) 


Note:  a  is  the  angle  between  the  0-x  axis  and  the  positive  tangential  direction 
of  coordinate,  s,  measured  counterclockwise 

Figure  1:  Sign  conventions  and  coordinate  systems 
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^  cosa  +  ^ILsina  +h  ^_  (1) 

9z  Bz  Bz  p  9z 

in  which  hp  is  the  distance  between  the  tangent  of  a  wall  element  to  the  centre  of  rotation 
’p’  of  the  cross  section  (it  can  be  proved  that  the  centre  of  rotation  is  also  the  shear  centre 
of  the  cross  section).  The  non-dependencies  of  a  and  hp  from  the  z  axis  in  this  equation 
indicate  that  the  beam  is  prismatic. 

The  equilibrium  equations  between  internal  resultants  and  external  forces  of  a 
cross  section  are: 

£  a  dAs  -  N  (a);  Jj  qhpds  -  Tp  (d) 
f  a  xdAc  -  M,  (b);  fqcosa  ds  -  Qt  (e)  (2) 

*/As  s  *  Js 

o  ydAs  -  My  (c);  Jqsina  ds  -  Qy  (0 

where  As  is  the  cross  section  area.  For  a  small  wall  element  as  shown  in  Fig.  1,  the 
equilibrium  equation  in  the  longitudinal  direction  of  the  beam  is 

*uiiLt-o  (3) 

9s  9z 


where  q  and  t  are  shear  flow  and  wall  thickness  respectively. 


The  compatibility  relation  between  displacement  D  along  the  s  axis  and  w  along 
the  z  axis  can  be  obtained  through  shear  strain  y  as: 


9w  9u 

_  +  _ 

9s  Bz 


(4) 


The  constitutive  relation  of  the  material  follows  Hooke’s  law  (linear  elastic 
material)  as 


7 


(5) 


where  x  and  G  are  sh„ar  stress  and  shear  modulus  respectively. 

2.2  Pure  Bending 

It  is  convenient  to  describe  the  geometry  of  a  cross  section  in  an  arbitrary 
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Cartesian  coordinate  system  (X,Y).  The  arbitrary  coordinate  system’s  origin  does  not 
necessarily  coincide  with  the  centre  of  gravity  of  the  cross  section  but  its  axes  are  parallel 
to  the  cross  section’s  x  and  y  axes  so  that  dX=dx=ds  cosa  and  dY=dy=ds  sina.  The 
transformation  between  the  two  Cartesian  coordinate  systems  can  be  obtained  through  the 
calculation  of  the  centre  of  gravity  as: 


x  -  X-Xc 


y  -  Y-Yc 


(6) 


By  assuming  no  shear  deformation  caused  by  transverse  shear  (Bemoulli-Naviers 
hypothesis),  using  equations  (1)  and  (3)  followed  by  integration  along  the  s  axis,  the 
longitudinal  displacement  w  can  be  expressed  as 

w  =  -Z,'X-t\'Y  +  wo(z) 

where  w0(z)  is  the  longitudinal  displacement  in  the  origin  of  tKe  system  (X,Y)  and  is  an 
unknown. 

Through  the  constitutive  relation,  the  axial  stress  caused  by  bending  can  be 
obtained  as: 

ab=  -E^//X-Er|//Y  +  Ew'(z)  (8) 

The  term  w0  can  be  eliminated  by  knowing  the  fact  that  there  is  no  axial  resultant  for  a 
beam  subjected  to  pure  bending.  Letting  the  axial  force  N  in  equation  (2a)  equal  zero  and 
using  the  coordinate  transformation  in  equation  (6),  yields: 

o  b  =  -E^'x-ErTy 

The  moments  of  inertia  of  a  plane  section  with  respect  to  the  x  and  y  axes  are  defined 
by  the  integrals: 


(10) 


The  moment  curvature  relationships  of  a  prismatic  beam  are  expressed  as: 

-E^  "j  x2dAs  -  Erj  "f  xydAs  -  Mx 
-Ei'jw.-En'iyVA.  -  My 
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By  substituting  these  quantities  into  equation  (9),  the  axial  stress  can  be  rewritten  as: 

(12) 


Ix-Iy  Iy-Ix 
a  „  -  JL _ +  _!1 _ 2_M 


D 


D 


where  D  =  1,1,  -  I  xy2. 


The  integration  of  equation  (2)  along  the  curvilinear  coordinate  leads  to 


$ 

qb(s,z)  =  qb*(z)- J 


9a 

9z 


tds 


(13) 


where  the  constant  qb‘  represents  the  shear  flow  at  the  beginning  of  the  integration.  The 
shear  flow,  qb*  vanishes  for  an  open  section  with  the  integration  starting  at  a  free  edge 
(shear  flow  equals  zero  at  this  point)  and  is  unknown  for  a  closed  section  in  which  the 
shear  flow  of  the  starting  integration  point  usually  is  non-zero. 

Introducing  two  integrals  which  represent  static  moments  of  a  portion  of  the  cross 
section  as: 


(14) 


and  substituting  equation  (12)  into  equation  (13),  one  can  obtain 


qb  “  qb- 


I  S  -I  S  I  S  -I  s 

y  *  *y  _  *  y  *y 


D 


■  n  *  y  *y 

D  y 


By  defining  two  quantities 


9b* 


D  ’ 


9by 


I  S  -I  s 

*  y  xy  » 

D 


(15) 


(16) 
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and  with  the  fact  that  qb  is  a  linear  combination  of  the  shear  forces  Q,  and  Qy,  equation 
(15)  can  be  written  as: 


I  S  -I  S  I  S -I  s 

(qb\-  ,..y  -  (qb*r  *  yp^y.  *Qy) 

(q:>q:)Q/(q;>q;)Qy 


(17) 


where  superscript  "u"  can  be  interpreted  as  shear  flow  caused  by  a  unit  of  shear  force. 
Since  there  is  no  twisting,  the  unit  angle  of  twist  is  zero  which  indicates  that 


Two  unknown  quantities  qbx*u  and  qby‘u  can  be  obtained  as 


They  are  equal  to  zero  for  an  open  section.  On  the  other  hand,  for  a  multi-cell  section 
with  n  compartments  these  equations  have  to  be  satisfied  for  each  cell  and  thus  result  in 
n  simultaneous  linear  equations. 

The  moments  of  inertia  in  equations  (12)  and  (15).  based  on  the  assumption  of  a 
thin  walled  section,  can  be  simplified  by  using  integration  by  parts  as: 

b  b 

Ix  -  J  x2dAs  =  Jx2tds  -  xSjo-Jsxdx  (20) 

1  o  o 


Because  the  first  term  on  the  right  hand  side  of  the  equation  vanishes  (according  to 
equation  (6),  x  is  normalized  so  Sx(0)  and  Sx(b)  equal  zero),  the  moments  of  inertia  can 
be  expressed  as 


2.3  Twisting 


There  are  no  lateral  displacements,  and  tj,  for  a  beam  subjected  only  to  pure 
twisting  around  its  centre  of  rotation.  Thus,  equation  (3)  can  be  reduced  to 


(22) 


According  to  the  basic  assumption  that  the  shear  deformation  of  the  cross  section  is  only 
caused  by  the  St-Venant  shear  stresses,  with  equation  (4),  yields 


dw 

"5T 


+ 


+ 


q5 


Gt 


(23) 


where  A  is  the  area  of  a  closed  section  surrounded  by  the  thin  walls.  The  second  term  in 
this  equation  vanishes  for  an  open  section. 

By  defining  dQ,  =  hpds  and  dQ2  =  J/(2At)ds,  integrating  equation  (23)  along  the 
s  axis,  gives 

w  =  -cp  'Qj  +  9  'Qj  +  wJz)  (24) 


The  warping  normal  stresses  can  be  obtained  by  using  Hooke’s  law  derived  as 

ou  -  -  Etp  "Q  j  +  Ecp  "Q2+Ew:(z)  (25) 

Letting  w07cp"=Q0,  and  following  Kollbrunner’s  [7]  definition  with  a  normalized  sectorial 
coordinate,  will  lead  to 

co  -  Q  .  -  Q  Q  (26) 

1  2  o 


allowing  the  axial  stress  in  equation  (25)  to  be  simplified  as 

o  -  -Ecp  "co  (27) 

The  quantity  Q0  in  equation  (26)  can  be  evaluated  using  the  fact  that  the  resultant  axial 
force  acting  on  the  cross  section  is  equal  to  zero.  Letting  N  in  equation  (2a)  equal  to  zero 
and  integrating  warping  normal  stresses  along  the  s  axis  yields 
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D 

J(Q  ,-Q  2)tds 


The  warping  shear  flow  can  be  obtained  from  equation  (3)  as 


s 

qw  -  Ecp  "'fa  t  dS  +  q* 


where  qu*  represents  the  shear  flow  at  the  beginning  of  the  integration  point. 

By  introducing  a  sectorial  static  moment  of  a  cut-off  portion  of  the  cross  section, 


to  tds 


the  warping  shear  flow  can  be  expressed  as 


-  E<P  "'S.tq: 

-  q<o°  +  qI 


The  warping  torsional  moment  with  respect  to  the  point  p  can  be  calculated 


T„  -  /qA*  '  Jq„^-dS-Jq„f2ids 

0  O 


The  assumption  (3)  that  the  longitudinal  shear  strain  is  only  caused  by  the  St-Venant 
shear  flow,  and  that  warping  shear  flow  has  no  contribution  to  the  shear  deformation, 
indicates  that 


/y.ds  -  J^Lds  -  0 
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Thus  the  last  term  in  equation  (32)  is  equal  to  zero 


(34) 


Using  integration  by  parts  of  the  first  term  in  equation  (32),  the  warping  torsional  moment 
can  be  obtained  as 


T.  -  <L<a 

b  b 

-  +  f_Lcotds  -  -Ecp  '"founds 

{  Zz  i 


or 

-  -E<p 

where  is  the  warping  constant 

b 

I.  “  Jw2tds 

o 


(35) 


(36) 


(37) 


1  I1V-  vv  OJ 


T 

_  0)  o 

%  ‘  ~j-S 


+  q<, 

/  OU  *  V,\*T' 

(qu  +  q^  )T» 


(38) 


The  constant  shear  flow  qu‘u  caused  by  a  unit  warping  moment  can  be  obtained  according 
to  assumption  (3)  or  equation  (33)  as 


q 


•  u 

(0 


(39) 
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Defining  M^T^,  the  warping  normal  stress  can  be  expressed  as 


The  derivation  of  the  above  equations  requires  that  the  centre  of  rotation  is  known 
in  advance  and  that  all  the  nodal  coordinates  refer  to  this  point.  Using  geometric 
relationships  leads  to 

h  ds  -  h  ds  -  x  sina  ds  +  y  cosa  ds  (41) 

pep  •'p 

Letting  d£2c  =  hc  ds  and  using  the  fact  that  dy=sina  ds  and  dx=cosa  ds,  the  normalized 
sectorial  coordinate,  co,  can  be  expressed  as 

0)  -  fl  -x  y+y  x-fl.-fl  (42) 

c  p-7  J  p  2  o 

The  centre  of  rotation  can  be  evaluated  by  using  the  fact  that  there  are  no  flexural 
moments  under  pure  twisting.  Letting  equations  (2b)  and  (2c)  equal  zero  gives 


b  b 

J o  xtds  -  0;  Jo  ytds  -  0 

O  0 

b  b 

or  Jco  xtds  -  0;  Jco  ytds  -  0 


Introducing  equation  (42)  into  (43),  gives 


b  b  b  b 

j  co  cxtds  -  xp  j  yxtds  +  ydJ*  x2tds  -  Jq  2xtds 

o  o  o  o 

b  b  b  b 

J  co  cytds  -  xp  Jy  2tds +  yp  f  xytds  -  2ytds 


Defining 


* 

Knc  -  JQ  cxtds;  hn,  “  Jfl2*tds 

o  o 

b  b 

c "  J^  cytds> 


ii 


and  introducing  these  quantities  into  equation  (44)  and  solving  two  simultaneous  linear 
equations,  the  centre  of  twisting  can  be  obtained  as 


yP 


^yO  c  ^yO  ^xO  c  ^xO  ^xy 

D 

^xO  c  "  Ko  ,)Iy  "  ^yO  c  "  ^yO  ,^xy 

D 


(46) 


The  quantities  in  equation  (45)  can  also  be  simplified 
as  equation  (20),  using  integration  by  parts  as 

b 

^xo  c  “  -Js.dQ  c; 

o 

k  "  "/Syd^  2;  " 


3.  Numerical  Procedure 

The  thin-walled  section  is  assumed  to  be  composed  of  a  number  of  narrow 
rectangular  segments.  Each  segment  is  numbered  consecutively,  starts  with  node  i  and 
ends  with  node  j.  The  coordinates  x,  y  and  0)  for  each  segment  are  distributed  linearly 
along  its  length;  therefore,  the  static  moments  S„,  Sy  and  S,,,  vary  parabolically  along  each 
element.  A  three  points  rule  of  Newton-Cotes  integration  can  be  used  to  obtain  an  exact 
solution  for  the  cross  sectional  constants  Ix,  Iy,  Ixy  and  I^.  The  integration  of  function  f(r) 
from  r=a  to  r=b  can  be  evaluated  as 

(f(r)dr  -  ^_i[f(a)+4f(il^)+f(b)]  <48> 

“62 


by  following  the  same  procedure 


•Js.dQ 

o 

>,dn: 


(47) 


With  the  equations  derived  in  the  previous  section,  the  step-by-step  procedure  to  evaluate 
the  cross  sectional  constants  is  as  follows: 

(1)  Input  nodal  coordinates  (X,,Yj)  and  (X^Yp  of  each  segment  k  with  respect  to  an 
arbitrary  coordinate  system  (X,Y),  its  thickness  t*  and  connectivity  with  other  elements. 

(2)  Calculate  coordinates  of  mid-point,  length  and  area  of  each  segment  k 

(3)  Obtain  centre  of  gravity  of  the  cross  section  as 


X  +x 

'  J 


\  - 

'  1  . 

2  * 

Y  .-t-  Y 

AXk  -  ) 

*k  - 

*  J  . 

-> 

AYk  -  ' 

Y-Y 
>  j 

Lk  -  J  AX2  +  AY2  ;  Ak  -  Lktk 


v  XXA.  v  IYA 

■'TaT  e’  Ta~ 

and  the  nodal  and  mid-point  coordinate  of  each  segment  (x,,y,),  (x^)  and  (xk,yk)  with 
respect  to  the  new  origin  (Xc,  Yc)  as 

x  -  X-Xc;  y  -  Y-Y 


(4)  Calculate  the  static  moments  Sx  and  Sy  of  a  cut  of  portion  of  the  cross  section  for 
points  i,  j  and  k  of  each  segment  as 


Ak(x(+  xk) 


4 

Vy.+y*) 


;  s. 


S^AS, 

Sy+ASy 


(5)  Obtain  the  moments  of  inertia  of  the  cross  section  as 


k  ax  *  AY 

ix  -  X^lsxW+4S^k)+sx(i)];  iv  -  4sy(iMSy(k)+sy(j)] 

n-l  6  n-1  O 

I„  -  t  4SA>-  s,(j)];  1„  -  t  s/i)- 4s,(k)- S,0)] 

O  n-l  o 


n- 1 


t6)  Calculate  AOck  as 


(7)  Evaluate  AQ^  as 


AOck  -  xkAYk-ykAXk 

AO  7k  -  2K-'Aii 
t„ 
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For  an  open  section,  this  term  is  equal  to  zero.  For  a  multi-cell  closed  section  with  n 
cells,  K  is  an  n  x  n  squared  matrix  as  described  in  Appendix  (A)  and  A  is  an  n  x  1 
matrix. 


(8)  Evaluate  Ix0c,  IyQc,  IxQ2  and  Iy£52  as 

I*.  -  Z^lS.OH4S,Ck)*S10)];  I 

o-l  O 

Cm  ’  Z^i[S,(i).4S,(kHS,(j)];  I,n, 

n-1  O 


*  AQ 

X— -[Sy(i)+4Sy(k)+S  (j)l 

D-  i  6 

*  AQ  lt 

X— -Ji[S  (i)+4Sy(k)+Sy(j)] 

n-1  6 


and  centre  of  rotation  of  the  cross  section  (Xp,yp)  as 


r  <I,n  c-  (I,QC- 

'  ‘  - D - ;  y>  ‘  " - D - 


(9)  Calculate  AQlk  as 


An,,  -  AQ  e,-x,AY,+  ypAX, 


(10)  Calculate 


A<n  -  An  An  2, 

(n,-n,),  -  (n^DjV.-Ain.-n,) 

Xl((n,*n!),*(n1.nI)i)^i) 

q  -  rJ _ _ 


(11)  Calculate  normalized  sectorial  coordinate  of  nodal  and  mid-point  of  each  segment 

to,  -  n,-fi2-oo 


(12)  Calculate  of  a  cut-off  portion  for  each  segment  as 

A. (co  .+ CD  .) 

AC  _  k  1  k  .  C  O 


;  S  -  S  +  AS 

’  CD  Cl)  Cl) 
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(13)  Obtain  the  warping  constant  1^  as 

I.  -  >  — — — -[S^^sjkys^ii)) 


'  AQ  ,-AQ  , 

Z - — 

U  6 


The  calculations  of  normal  stresses  are  straight  forward  and  do  not  need  to  be 
discussed.  Only  the  procedure  for  calculating  the  transverse  and  the  warping  shear  stresses 
are  presented  here: 


(14)  Calculate  q^0"  and  q^"1  as 


qb* 


- - - , 


Qby  ” 


D 


(15)  Calculate  the  total  shear  flow  around  each  cell  as 


(O.  -  (Z— iq«“"(i)-4q17<k).qir(j)]]..ll 


«><„,  -  G—  !<<■)»  4q~(k)  H- 

lk 


cell 


and  the  two  quantities  qb/u  and  qbv*u  as 

<q*>  -  -[Kr'kqDJ 
v  -  -iK]-H(q“)  j 

(16)  The  transverse  shear  flow  caused  by  a  unit  shear  force  can  be  obtained  as 

„  u  ou  »u 

qb»  “  qbx  +  qbx 

_  u  ou  «u 

qby  “  qby +  qb> 

(17)  Calculate  qj"3  as 


OU 

qu 
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(18)  Calculate  the  total  shear  flow  around  each  cell  as 


and  qM’u  as 


/  OUN 

(q« ) 


cell 


'(i)  +  4q:u(k)  +  q;u(j)]] 


cell 


W  -  'IKIIOJ 


(19)  The  warping  shear  flow  caused  by  a  unit  shear  force  can  be  obtained  as 

u  ou  _  *  u 

q<u  —  q<i) +  9o) 

(19)  The  shear  flows  caused  by  different  forces  now  can  be  evaluated  by  multiplying 
these  unit  shear  flows  with  applied  forces. 

The  above  procedure  can  be  applied  on  not  only  the  closed  sections  but  also  the  open 
sections  by  letting  the  quantities  Q2  and  9*  equal  zero. 
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4.  Examples 

Two  examples  are  presented  in  this  section  in  tabular  form  in  order  to  demonstrate 
the  above  procedure.  The  Arabic  numeral  in  the  first  row  of  each  table  indicates  the  step 
number.  The  results  can  be  verified  with  the  equilibrium  condition. 

Example  1:  A  single-cell  box  beam  with  different  wall  thicknesses  is  shown  in  figure 
2.  (a)  Calculate  the  cross  sectional  constants  and  shear  centre,  (b)  Calculate  the  transverse 
shear  flow  caused  by  transverse  shear  forces  Q,=-5.5  x  10s  kN  and  Qy=-4.125  x  105  kN. 
(c)  Calculate  the  warping  shear  flow  due  to  a  unit  warping  moment. 


Transverse  shear  flow 


Warping  shear  flow 


(Normalized) 

Figure  2:  Single-cell  box  beam 


17 


21  18.6 


-341.06  -95.75 


10004  94.32 


88.6 


-D6  56.84  15.91 


-16.67  -15.72 


(2) 


(3) 


(4) 


(5) 


(8) 


ASy  Sy  I  1+4+1  I  Iy  I  Iy*  I  lyffc  I  Iy 


105  I  1 07  10’  10’°  !0'° 


elem  ty  node  Y  AY  Ly  Ak  YkAy 

105  10s 


0 


1  10 _  150  300  500  5  7.5  -148.6  _ |  -5.59 1 -29.79 1  -89.37 1 -1 19. 16f  -27.66|  -33.1 

-1.84 


2 

300 

r 

1.4 

420  240  300  3  12.6  121.4  -6.51  -37.26  -89.42  67.07  -27.94-24.84 


2.72 


241.4  -3.79 


8.32 


390  -300  500  10  39  91.4  4.53  19.68  -59  04  -78.72  11.24  10.93 


0.82 


240 

-58.6 

I 


120|  -240 1  300 1  3 1  3.6|  -178.6  3  57 1  19.63)  -47.1 1 1  35.33|  14  72|  13.09 


-3.56 


0 


0 

-298.6 

z 

21  62.7 

-284.94 

-95.48  -29.64  -33.92 

298  6 

-E/6  47.49 

15.91  4.94  5.65 

19 


21 


Example  2:  For  a  two-cell  section  with  constant  wall  thicknesses  as  shown  figure  3, 
calculate  (a)  the  cross  sectional  constants  and  shear  centre,  (b)  the  transverse  shear  flow 
caused  by  Q,  and  Qy  equal  to  105  kN,  and  (c)  the  warping  shear  flow  caused  by  a  unit 
warping  moment. 


[K] 


216.56  -56.56 
-56.56  216.56 


(A)  - 


I,  =  287.66  x  107 
Iy  =  287.66  x  107 
Ixy  =  -53.69  x  107 
I<o  =  2453.36  x  1010 
J  =  2.88  x  109 


24x1 04 
24x1 04 


w  -  (q;u) 


-2987x10  s 
15887xl0'5 


Xc  =  327.5 
Xp  =  -65.04 
Xp  =  262.46 


Yc  =  327.5 
yp  =  -65.04 
Yp  =  262.46 


5217.33xl0"5 

52l7.33xlO'5 


Figure  3:  Multi-cell  section 
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5.  Concluding  Remarks 

This  report  has  presented  the  complete  mathematical  derivations  of  the  equations 
used  for  the  calculation  of  cross  sectional  constants  and  stress  distributions  of  thin-walled 
sections.  The  derivations  are  based  on  the  same  assumptions  as  Kollbrunner’s  [7]  open 
section  and  extend  to  multi-cell  closed  sections.  These  equations  can  be  used  for  any 
shape  of  cross  section,  both  open  and  closed,  and  are  independent  of  the  orientation  of 
the  cross  section.  A  numerical  procedure  based  on  these  equations  has  been  outlined  and 
verified  with  several  examples.  It  can  be  easily  implemented  into  the  current  SCRAP 
program. 

It  should  be  noted  that  these  equations,  which  are  based  on  certain  assumptions, 
have  some  limitations.  The  major  two  are  that  the  derivations  deal  only  with  prismatic 
members,  and  cross  sections  retain  their  shape  during  deformation.  Thus,  the  solutions  of 
the  equivalent  beam  models  can  only  be  improved  with  increasing  numbers  of  beam 
elements,  and  are  only  valid  for  the  unbuckled  state  of  the  ship  hull’s  elastic  response. 
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Appendix  A  -  St-Venant  Torsion 

The  governing  differential  equation  of  a  prismatic  beam  subjected  to  St-Venant 
torsion  is 

T  -  GJtp  '  (Al) 


where  J  is  torsional  constant. 

For  an  open  thin-walled  section  which  is  composed  of  a  number  of  thin 
rectangular  elements,  the  torsional  constant  can  be  computed  as  the  sum  of  the  values  for 
the  individual  element. 


j  -  yibt  (A2) 

_ .  ^11 

where  i  is  the  i^  element.  The  shear  stress  is  linearly  distributed  across  the  thickness  of 
the  wall  w'ith  a  zero  average.  Its  maximum  value  of  shear  stress  is  at  the  wall  surface  and 
can  be  written  as 


Tt 

(x  ,),  -  -11 

J 


For  a  closed  thin-walled  section  with  multi-cells,  the  torsional  constant  is 

J  -  4(A)t[K]'1(A} 


(A3) 


(A4) 


where  A,  is  the  area  of  cell  i  and  the  entries  in  matrix  [K]  are 


k 


ii 


ds 

IT  T 


(A5) 


The  subscript  "ii"  represents  a  summation  performed  along  walls  of  cell  1  and  ”ij" 
represents  a  summation  performed  along  common  w'alls  between  cell  i  and  j. 

For  a  single  cell  section,  equation  (A4)  can  be  reduced  to 


J 


4A2 


(A6) 
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For  a  beam  subjected  to  a  specified  torsion  T„  the  rate  of  twisting  angle  <p'  can  be 
obtained  from  equation  (Al).  The  St-Venant  shear  flow  (qj,  in  each  cell  then  can  be 
found  as 


(qs)  -  2Gtp  '[Kj-'lA) 
-  IllKV'iA) 


(A7) 
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Appendix  B  -  Open  Section  -  Example 


For  a  channel  section  with  constant  wall  thickness  as  shown  below,  calculate  the 
shear  centre  and  warping  constant  and  verify  the  results. 


The  exact  solutions  of  the  shear  centre  and  warping  constant  of  a  channel  section  with 
constant  wall  thickness  are 
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Xc  =  151;  Yc  =  215 
xp  =  118.85;  yp  =  -160.16 
Xp  =  269.85;  Yp  =  54.85 
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